ABSTRACT
Background: Coronavirus disease 2019 (COVID-19) is a disease that causes fatal disorders including severe pneumonia. Our study aimed to utilize bioinformatics method to analyze the expression profiling by high throughput sequencing in human bronchial organoids/primary human airway epithelial infected with SARS-CoV-2 to identify the potentially crucial genes and pathways associated with COVID-19.Methods: We analyzed microarray datasets GSE153970 and GSE150819 derived from the GEO database. Firstly, the Differentially expressed genes （DEGs） in human bronchial organoids/primary human airway epithelial infected with SARS-CoV-2. Next, the DEGs were used for GO and KEGG pathway enrichment analysis. Then, the PPI network was constructed and Cytoscape was used to find the key genes.Results: Gene expression profiles of GSE153970 and GSE150819, in all 12 samples were analyzed. A total of 145 DEGs and 5 hub genes were identified in SARS-CoV-2. Meanwhile, we found that the 145 genes are associated with immune responses and the top 5 hub genes including CXCL8, CXCL1, CXCL2, CCL20, and CSF2 were mainly related to leukocyte migration, endoplasmic reticulum lumen, receptor ligand activity. In addition, the results also showed that the hub genes were associated with Cytokine−cytokine receptor interaction, IL−17 signaling pathway, and Rheumatoid arthritis in SARS-CoV-2 infection.Conclusion: The five crucial genes consisting of CXCL8, CXCL1, CXCL2, CCL20, and CSF2 were considered as hub genes of SARS-CoV-2, which may be used as diagnostic biomarkers or molecular targets for the treatment of SARS-CoV-2. It is evidenced that bioinformatics analyses in SARS-CoV-2 can be useful for understanding the underlying molecular mechanism and exploring effective therapeutic targets.